Scaling behaviour in the fracture of fibrous materials 
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We study the existence of distinct failure regimes in a model for fracture in fibrous materials. 
We simulate a bundle of parallel fibers under uniaxial static load and observe two different failure 
regimes: a catastrophic and a slowly shredding. In the catastrophic regime the initial deformation 
produces a crack which percolates through the bundle. In the slowly shredding regime the initial 
deformations will produce small cracks which gradually weaken the bundle. The boundary between 
the catastrophic and the shredding regimes is studied by means of percolation theory and of finite- 
size scaling theory. In this boundary, the percolation density p scales with the system size L, which 
implies the existence of a second-order phase transition with the same critical exponents as those of 
usual percolation. 
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I. INTRODUCTION 



The nature of fracture of non-homogeneous materials 
is an important problem in material science research. 
Computer simulation of the fracture phenomenon are 
very useful since the analytical approachs is very difficult 
to perform. This difficulty arises from the non-uniform 
character of the material and their discrete nature, which 
are fundamental ingredients for understanding the rup- 
ture process Usually, computer simulation in these 
materials gives interesting results, however the high de- 
gree of correlations between the constituents leads to a 
high computational cost. Bundles of unidirectional fibers 
form a system with low degree of correlations allowing the 
fracture process be simulated in a large scale. 

The study of fibrous materials is not recent, as one can 
find in the work of Daniels [^, who in 1941 studied the 
rupture of a bundle of fibers with a known probability 
distribution of strength. Recently, Hansen and Hemmer 
1^ studied the distribution H{S) of the sizes S of burst 
avalanches, i. e., an instantaneous propagation of a crack. 
They found a power-law behaviour: H(S) ^ S~°' with 
the exponent a depending on how the load is shared be- 
tween the fibers. For global load sharing, where the load 
is shared equally among non breaking fibers, they ob- 
tained a = 2.5. For local load sharing, when a fiber 
breaks its load is shared among nearest-neighbours no 
breaking fibers, they obtained a = 4.5. Those results 
have been obtained for a one-dimensional lattice of fibers, 
that is, the load sharing is the only correlation between 
the fibers. Other approaches have been introduced to 
discuss this problem MM- The role of the homogeneous 



support matrix on the failure of composite materials has 
also been discussed by some authors |^-|^. 

A model of a bundle of unidirectional fibers, which 
takes into account external parameters like temperature 
and velocity of traction, has been proposed in 1994 by 
Bernardes and Moreira In this model, the correla- 
tions between the fibers are present through the proba- 
bility of rupture of a fiber, which depends on the num- 
ber of unbroken nearest neighbouring fibers. A cascade 
mechanism - inspired on models for avalanches - is used 
to propagate cracks through the material: when a fiber 
breaks, its neighbours are visited and can break too. In 
this model, all unbroken fibers have the same deforma- 
tion, i.e, one has global load sharing. However, the cas- 
cade mechanism introduces a local effect. In a subsequent 



work 1 10 1, the dependence of the frequency of cracks with 
the crack sizes were used to determine the failure regimes. 
Two basic regimes were discussed: a regime where cracks 
of the size of the system were present and another one 
where only small cracks appeared. Those regimes were 
identified, respectively, with the brittle and ductile fail- 
ure regimes. The criterion used to distinguish one regime 
from another was based on self-organized criticality, i.e, 
in the brittle-ductile transition region, cracks of all sizes 
were present. However, they did not take into account 
finite size effects in their analysis, which are very impor- 
tant in this type of process. In fact, a finite-size scaling 
analysis should be performed, in order to guarantee a 
better definition of the failure regimes. 

The aim of the present paper is to introduce a criterion 
which defines the failure regimes on fibrous materials. 
When a fiber bundle breaks, two regimes can be observed: 
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a catastrophic regime, when a sufhciently large number 
of fibers are simultaneously broken, and a slow process of 
successive rupture of fibers, here called shredding regime. 
The first regime occurs at low temperatures and/or high 
strains and are similar to a brittle fracture. It is charac- 
terized by the fact that an initial deformation produces 
a large crack which percolates through the bundle. The 
shredding regime occurs for higher temperatures and/or 
lower strains, and is similar to the ductile regime. In this 
case, the first deformations produce small cracks which 
weaken the bundle and thus cause its failure. The crite- 
rion is implemented by considering the static failure of a 
modified version of the model introduced by Bernardes 
and Moreira . A second order phase boundary between 
two regimes is found for a given strain. A finite-size scal- 
ing analysis is used to determine the critical temperature 
and exponents. 



II. THE MODEL 

The model for the fibrous material here discussed con- 
sists of a bundle of iVo = L x L parallel fibers with a 
cross-section forming a triangular lattice. Each fiber has 
the same elastic constant fc, and they are fixed at both 
ends to parallel plates. One plate is fixed and the other 
plate can be pulled by an external force. When the bun- 
dle is pulled by a force F, all fibers undergo the same 
linear deformation z — F/Nk, where N is the number 
of unbroken fibers. We assume that a fiber has a fail- 
ure probability which increases with the deformation z. 
When this deformation reaches a critical value z^, the 
breaking probability of an isolated fiber is equal to one. 
When the bundle has a deformation z, a fiber i has a 
failure probability related to its elastic energy and to the 
number of unbroken neighbouring fibers n^, given by 



On = — = 



Fn 



Z/Zc 



1 



exp 



(fczV2) - (fczg/2) 
KbT 



(1) 



Defining the strain of the material as 5 — z/ Zc and the 
normalized temperature as t = KbT/E^, where T is the 
absolute temperature, Ec = is the critical elastic 

energy and Kb is the Boltzman constant, we can rewrite 
the failure probability as 



■ exp 



(2) 



This definition of the failure probability is different 
from that used by Bernardes and Moreira |9[], since now 
we have introduced (5 as a multiplicative factor to impose 
that, for (5 = 0, P^{S) = 0. 

The static failure of a fiber bundle is produced by ap- 
plying a constant force Fq to the bundle, for example, by 
hanging a weight on the moving plate. The initial strain 
of the bundle is given by 



Nokzc 



(3) 



The simulation of the rupture process proceeds as fol- 
lows. At each time step of the simulation, we randomly 
choose a set of Nq{= qNo) unbroken fibers, where the 
number q represents a percentage of fibers and it allows 
us to work with any system size. So, differently of an 
Ising model, where all the sites are "tested" at each time 
step, in our model only a number Nq of randomly chosen 
unbroken fibers are tested. It represents the continuous 
growth of the bundle due to the continuous traction. For 
each chosen fiber, we evaluate the probability of rupture, 
using Eq. 1, and compare it with a random number in 
the interval [0,1). If the random number is less than 
the failure probability, the fiber breaks. To simulate the 
load spreading, the same process is repeated for all neigh- 
bouring unbroken fibers. The failure probability of these 
neighbouring fibers increases due to the decreasing of rii 
and a cascade of breaking fibers may begin. This pro- 
cedure describes the propagation of a crack through the 
fiber bundle, which occurs in all directions perpendicular 
to the force applied to the system. The cascade process 
stops when the test of the probability does not allow the 
rupture of any other fiber on the border of the crack or 
when the crack meets another already formed crack. This 
collision leads to the fusion of cracks, and it is the mech- 
anism to explain the rupture of the material in the shred- 
ding regime. The same cascade propagation is attempted 
by choosing another fiber of the set A'',. After all the A^^^ 
fibers have been tested, the strain is increased if some 
fibers have been broken. This new strain is the same for 
all the remaining unbroken fibers. Since the force is fixed 
(the weight hung on the bundle), the greater the num- 
ber of broken fibers, the larger is the strain on the fibers, 
and the higher is the failure probability. Then, other 
set of Ng unbroken fibers are chosen and the rupture 
process restarts. This process stops when all the fibers 
are broken, i.e, the bundle breaks apart. In this model, 
a combination of local and global load sharing occurs. 
That is, after a fiber breaks, a cascade may begin which 
simulates the local load sharing. When the cascade pro- 
cess stops, the stress is distributed equally between all 
unbroken fibers which is the global load sharing. 



III. RESULTS 

The failure probability (Eq. |^) can be written as 



yrii + 1) 

where we introduce the parameter r(t, 5) defined as 



T{t,S) = 5 exp 



t 



(4) 



(5) 
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For a triangular lattice (with coordination number 6) 
and T{t, 5) > 6, the rupture of any fiber induces the rup- 
ture of the whole bundle, i.e., the bundle breaks with just 
one crack. Obviously, this crack forms a cluster which 
percolates through the entire system. 

We can define the density of the percolation crack as 



where Npc is the number of broken fibers belonging to 
the percolating crack. Thus, when T{t, 5) > 6, we have 
p = 1. On the other hand, as it has been observed in 
previous works pO| , pl| , for higher temperatures and/or 
lower strains, the fracture of the bundle is caused by 
many small cracks, none of then large enough to perco- 
late through the system. Thus, for a fixed temperature, 
if one starts with a large enough strain do and one de- 
creases it, the system goes from a regime where p = 1 to 
another regime where p — » 0. This behaviour is the same 
as the one encountered in the percolation problem. 

Figure 1 shows the density of the percolation cluster p 
versus the initial strain Sq, for two different temperatures. 
As one sees, p = 1 for high values of Sq, and jumps to 
zero for low enough value of ^o- So, we may assume that, 
for a fixed temperature, there is a critical value Sqc above 
which one observes a percolation crack, and below which 
there is no percolation at all. Another interesting feature 
that one can observe in Figure 1 is that, if one substitutes 
into Eq. ^ the values of t and Sqc corresponding to the 
transition region {Sq ^ 1.18 for t = 1.0 and 6o ^ 1.37 
for t = 4.0), we get for both instances T{t,d) ^ 1.73. 
The fundamental reason for obtaining this value wil be 
explained below. 

In contrast to that described above, the same be- 
haviour does not occur when we keep 6o fixed and change 
the temperature. Figure 2 shows the results obtained for 
the density of the percolation cluster p versus tempera- 
ture t, for ^0 = 1-4. We observe that, initially, p decreases 
as the temperature t increases, and around t ^ 4.5, the 
value of p seems to go to zero. However, an additional 
increase in the temperature will revert the process and a 
minimum appears. Note that at the point of minimum 
again (r{t = 4.5, i5 = 1.4) ~ 1.73), which is the same 
value reported above. For low temperatures {t < 2.0), 
when a fiber breaks, the probability is so high that this 
rupture initiates a cascade which breaks the whole bun- 
dle. By increasing the temperature, a number of small 
cracks are formed, inhibiting the formation of a perco- 
lating cluster and the density p decreases. However, all 
those processes occur in the first step of the simulation 
when Nq attempts to break the bundle are performed. 
Thus, for t < 4.5, the bundle has been broken due to 
the crack which percolates the system during the first 
Nq attempts to break it. For t > 4.5, all the first Nq 
attempts do not succeed to generate a crack which per- 
colates the bundle. However, some fibers have been bro- 
ken and cracks were formed. In the second step of the 



simulation, a new value of 6 is used (higher than Sq) and 
a new set of Nq trials are chosen. But now one has a 
higher value for T{t, 5) therefore it is easier to produce a 
large crack which percolates the bundle. By increasing 
the temperature, a smaller number of fibers arc broken 
in the first Nq attempts, and then, in the second step, 
the are more unbroken fibers and therefore the density p 
increases, thus forming a minimum in the graph of Figure 
2. 

In fact, we can assume that there is a critical value 
for T{t, S) ~ 1.73 that defines the transition between two 
regimes. In the first one, a catastrophic fracture occurs 
due to the first attempt to break the bundle, while in 
the second case the rupture of the bundle occurs due to 
the formation of small cracks, which weaken the bundle. 
A percolating crack may also occur in the second case, 
however the fracture dynamics is given by the weakening 
of the bundle not by the catastrophic propagation of a 
crack. 

In order to consider the present model in the context 
of percolation theory, we shall use the parameter F as an 
arbitrary parameter without regarding it as a function of 
the strain S and temperature t. Within the percolation 
point of view, we map the original model into a triangular 
lattice where the empty sites corresponds to the unbroken 
fibers. The parameter F is the analog to the percolation 
probability. The algorithm for the mechanism of fracture 
is mapped into the following algorithm for the percola- 
tion problem. An empty site (unbroken fiber) i is chosen 
at random; Its occupation (failure) probability Pi is cal- 
culated by dividing F by the number of its neighbouring 
empty sites (unbroken fibers) plus one. This probabil- 
ity Pi is compared with a random number r G [0, 1); If 
Pi > r, the site is occupied (the fiber is broken) and a 
cluster (crack) can be formed, i.e, an empty neighbouring 
site (an unbroken neighbouring fiber) is randomly chosen 
and the process are repeated; Otherwise, another site, on 
Nq in total, is chosen. 

When a cluster is formed, we test if it percolated 
through the system. If it does, we calculate the density p 
of the percolating cluster. Figure 3a shows the results ob- 
tained for several system sizes. Two regions are separated 
by the transition point Tc- The larger the system size, 
more clearly is the transition between those two regions. 
Observe in the detail, shown in Figure 3b, that a second 
order phase transition takes place at Tc = 1.733(1). This 
implies that, at that point, clusters of all sizes should be 
present, as confirmed by the results shown in Figure 4. 
In this figure, the results have been obtained for a system 
size L = 5000 (2.5 x 10^ fibers) and averaged over 1000 
samples (it took nearly 24h on a Sun Enterprise 8GB 
computer) which gives the following power law 

H{S)^S-^ , (7) 

where r — 2.037±0.007. A finite size analysis can be per- 
formed by plotting t{L) as a function of L~^^'^, where v 
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is the exponent related to the divergence of the correla- 
tion length at the transition. We tested several value of v 
and the best linear fitting were obtained for — 4/3, as 
shown in Figure 5. This value corresponds to the exact 
exponent v for percolation at d — 2. The value of the 
exponent r for an infinite lattice is, then, evaluated to be 
t(oo) — 2.05 ± 0.01, in an excellent agreement with the 
theoretical value, Toq = 2.055 

In order to check if our problem belongs to the same 
universality as the percolation problem, we have done a 
finite-size scaling analysis by assuming the scaling law 



where 



(8) 



(9) 



is a universal function of eL^^'^ only, and /3 and v 
are the critical exponents for the infinite lattice. Figure 
6 shows the finize-size scaling plot pL^^'^ versus eL^^" for 
nine sizes of L. We have used Fc = 1.733, i' — 4/3 and 
the best value of /? which validates Eq. || is /? = 0.14. 
This value is also in an excellent agreement with the 
known value for the usual percolation. 

Now, returning to the original fracture model, we use 
Eq. |5| to obtain the critical temperature tc in terms of 
the critical parameter F^ and of the initial strain Sg 



tr. = 



51 



1 



In(Fc) - ln(5o) 



Using this expression we plot the fracture regimes di- 
agram, in the temperature t versus the initial strain (5o 
plane, depicted in Figure 7. Two fracture regimes are 
separated by a second order transition line. In region 
C the fracture is catastrophic and in region S we have 
the shredding regime. Note that the catastrophic regime 
only occurs for (5o > 1 and for low temperatures. In this 
figure, the solid line corresponds to the analytical results, 
and the points were obtained by simulations. 



IV. CONCLUSIONS 

In conclusion, we have studied a model for fracture in 
fibrous materials in (2-|-l)-dimensions and shown the ex- 
istence of two failure regimes: the catastrophic regime, 
where the initial deformation produces a single crack 
which percolates through the bundle; and the slowly 
shredding regime, where the initial deformation produces 
small cracks which gradually weaken the bundle. By us- 
ing percolation theory and finize-size scaling arguments, 
we were able of finding the transition line between these 
regimes. Our results indicate that this transition is of 



second order. Finally, we have shown that this model 
belongs to the same universality class as the percolation 
problem. 

Acknowledgements We thank H. J. Herrmann and J. 
A. Plascak for fruitful discussions and suggestions. C. 
Moukarzel and M. Continentino did important sugges- 
tions about the approach through percolation theory. We 
also thank O. F. de Alcantara Bonfim for helpful criti- 
cism of the manuscript. One of us (ATB) acknowledges 
the kind hospitality of the Depto de Fi'sica, UFMG. We 
also acknowledge CNPq and FAPEMIG (Brazilian agen- 
cies) for financial support. Most of our simulations have 
been performed on the Sun Enterprise 8GB computer of 
the CENAPAD MG/CO. 



[1] 



(10) [6; 



[7] 
[8] 

[9: 

[lo; 
[11 

[12 



Statistical Models for the Fracture of Disordered Media, 
edited by H. J. Herrmann and S. Roux, (North-Holland, 
Amsterdam, 1990). 

Daniels, H. E., Proc. R. Soc. A 183, 404 (1945); 
Hansen, A. and Hemmer, P.C., Phys. Lett. A 184, 394 
(1994) 

Duxbury, P. M. and Leath, P. L., Phys. Rev. B 49, 12676 
(1994); Leath, P. L. and Duxbury, P. M., Phys. Rev. B 
49, 14905 (1994) 

Zhang, S. D. and Ding, E. J., Phys. Lett. A 193, 425 
(1994); J. Phys. A 28, 4323 (1995); Phys. Rev. B 53, 
646 (1996) 

Ibnabdeljalil, M. and Phoenix, S. L., Acta Metall. Mater. 
43, 2975 (1995); Beyerlein, I. J., Phoenix, S. L. and Sas- 
try, A. M., Int. J. Solid Structure 33, 2543 (1996) 
Zhou, S. J. and Curtin, W. A., Acta Metall. Mater. 43, 
3093 (1995) 

Moraes, W. A., Godefroid, L. and Bernardes, A. T., Proc. 
of the Sixth International Fatigue Congress. Berlin, vol 
HI, 1499 (1996) 
Bernardes, A. T. and Moreira, 
15035 (1994) 

Bernardes, A. T. and Moreira. 
5, 1135 (1995) 

Menezes-Sobrinho, I. L., Moreira, J. G. and Bernardes 
A. T., Int. J. Mod. Phys. C 9, 851 (1998) 
StaufFer, D. and Aharony, A., Introduction to Percolation 
Theory (Taylor & Francis, London, 1992) 



J. G. 



J. G. 



Phys. Rev. B 49, 
J. Phys. I (Paris) 



February 1, 2008 



4 



IV CONCLUSIONS 





1.70 



<1 <iL = <IUIII 



1.75 

r 



1.80 




FIG. 1. Density of the percolating cluster p versus the ini- 
tial strain So for two different temperatures: t = 1.0 (filled 
circles) and t = 4.0 (open diamonds). The system size is 
L = 1000 and the data were averaged over 1000 statistically 
independent samples. 



FIG. 3. (a) Density of the percolating cluster p versus F 
for five different system sizes; (b) Zoom of the region cor- 
responding to the small box drawn in the left plot, showing 
more clearly the crossing of the curves. The data were aver- 
aged over 1000 statistically independent samples. 




FIG. 2. Density of the percolating cluster p versus temper- 
ature for an initial strain So = 1.4 and three different system 
sizes: L = 800 (open squares); L — 900 (filled circles) and 
L = 1000 (open diamonds). The data were averaged over 
1000 statistically independent samples. 




FIG. 4. Log-log plot of the averaged number of cracks H{S) 
versus the crack size 5 for L = 5000 at Fc = 1.733. The points 

have been obtained by averaging over 1000 statistically in- 
dependent samples. The data show a power law behaviour 
(expected at the criticality) with exponent r — 2.037 ±0.007. 
The insert shows a detail of the whole set. The solid line in 
this insert has exponent 2.037. 
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FIG. 7. Fracture regimes diagram of the temperature t in 
function of the initial strain 60, where C represents the catas- 
trophic regime and S represents the shredding regime. Sohd 
hne represents the theoretical curve and filled circles represent 
the data obtained in our simulations. 
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